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Abstract. The problem of electron scattering on the one-dimensional complexes is considered. We propose 
a novel theoretical approach to solution of the transport problem for a quantum graph. In the frame of the 
developed approach the solution of the transport problem is equivalent to the solution of a linear system of 
equations for the vertex amplitudes 'S'. All major properties, such as transmission and reflection amplitudes, 
wave function on the graph, probability current, are expressed in terms of one F-matrix that determines 
the transport through the graph. The transmission resonances are analyzed in detail and comparative 
analysis with known results is carried out. 



PACS. 73.23.-b Electronic transport in mesoscopic systems 
materials and structures - 85.35.Gv Single electron devices - 



- 73.63.-b Electronic transport in nanoscale 
03.65.Nk Scattering theory 



1 Introduction 

Impressive progress in the development, manufacturing 
and understanding of various nanostructures is observed 
in recent times [Tl[21llll[51[Sl[71[Hl[SJ[in] . This makes scien- 
tists not only to formulate the new theoretical problems 
and to reconsider the existing ones but also to look for the 
new approaches to the solution of these problems. That 
"new" approaches usually have their own history. 

The classical example is the dynamics of the electron in 
a graph (one-dimensional complex JT^). This problem has 
more than seventy years history and was introduced into 
physics in connection with calculation of the properties 
of conjugated organic molecules P^[T^[T^I151I161[T7] . but 
then gradually lost its relevance to the quantum chemistry. 
Nevertheless, the model of electron in a graph stimulated 
the emergence of the new field in the mathematics related 
to the investigation of the spectrum of self- adjoint opera- 
tors in compact and noncompact graphs 18,19,20,21,22, 
123) . In recent years the renewed interest in this area in- 
volved significant redevelopment and enhancement of the 
existing theory [Mll^l^ . 

The quantum mechanics in the graph demonstrates 
many interesting peculiarities. One of them is the spec- 
trum of the Schrodinger operator. To find it one should 
construct its self-adjoint extension which is not unique |19| . 
This is most commonly done in way that was first pro- 
posed by Griffith [37]. The classical counterpart of the 
Griffith's extension is the well known Kirchhoff rules for 
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current conservation in each node supplemented with the 
wave function continuity condition resulting from the quan- 
tum mechanics ideology. The spectrum of this self-adjoint 
extension is discrete on compact graph [221128] . 

The more interesting property of the spectrum is the 
existence of the discrete eigenenergies embedded in the 
continuum (in what follows, it is abbreviated as BIG). 
The possibility of such states was shown in the promi- 
nent paper of Wigner and von Neumann [29| and the 
method to construct them was proposed. Later, it was 
improved by Stillinger and Herric [30]. The central idea 
of this approach consists in such a choice of the potential 
that the continuum wave function becomes integrable. In 
the same work of Stillinger and Herrick the BIG were pre- 
dicted in the semiconductor super lattices. In 1992 it was 
confirmed in the experiment on the IR irradiation of the 
Alo.48lno.52As/Gao.47Ino.53As snperlattice [31j . 

The discrete spectrum superimposed on the continuum 
in quantum graph was considered in details in a number of 
papers [Ml[5^[55] . Such states in the graph behave quite 
similarly to the ordinary (in the presence of potential) 
BIG. This fact encouraged Bulgakov et.al. 32 to use the 
term BIG in the case of graph. In the works of Texier [211 
155] the presence of the discrete eigenvalues for positive en- 
ergies on the graph is more carefully interpreted as decou- 
pling of the compact part of the graph from the external 
leads. 

The considered phenomena can be interpreted in the 
formalism of the de Broglie waves dynamics [34j . However 
the results of Xia [33] are rather cumbersome to analyze 
and to make any generalizations. 
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In connection with the above the aim of the presented 
work consists in the development of the successive mathe- 
matical technique based on the scattering theory to solve 
the following problems: wave function of the electron in 
the one-dimensional complexes (graph-like structures) , en- 
ergy spectrum, transmission and reflection amplitudes, 
Green's function and response to the external electromag- 
netic field. These objects are treated in the frame of the 
single technique in contrast to the current papers in this 
field that use the different mathematical methods. In addi- 
tion, the developing approach gives the results with clear 
physical meaning. Moreover, using the proposed technique 
allows not only to explain the some fascinating effects (24l 
[?T1[5^ , but also to formulate the conditions of their occur- 
rence. 

This paper is organized as follows. In section [2] we 
briefly consider the basics of electron transmission through 
the quantum graphs; we give the general considerations 
and equations for electron moving on the quantum graph. 
The main mathematical formalism of the developed method 
and its application to the electron transmission through 
the potential with finite support is described in section [31 
In section|3]we apply the proposed formalism to the prob- 
lem of electron transmission through symmetric and asym- 
metric (with and without Aharonov-Bohm flux) quantum 
ring. In section [S] the obtained results are summarized. 
The application of the vertex amplitudes method to the 
classical quantum mechanics problem of electron scatter- 
ing by the potential with finite support is presented in 
Appendix 1X1 In Appendix [Bl the problem of electron scat- 
tering by the system of parallel coupled quantum wells is 
considered. 



2 Preliminaries 

Most of the problems of one-dimensional quantum me- 
chanics can be formulated in the language of quantum 
graphs. Moreover, the solution of these problems in the 
frame of the graphs formalism could be more efficient 
and "elegant" in comparison with the traditional meth- 
ods of quantum mechanics. One of such problems is the 
electron scattering on the system of connecteed in paral- 
lel quantum wells (the detailed analysis is carried out in 
Appendix [Bj . A simple and clear example which demon- 
strates the use of the graphs formalism is the problem of 
the one-dimensional quantum well (see fig.[l]). 

There are various ways to solve the Schrodinger equa- 
tion on the graph |34II241[^ . In the presented work we 
develop the vertex amplitudes method. The main idea of 
this method is to express the parameters of the problem 
through the values of the wave function at the vertex of 
graph "ifi and '52. As will be shown below the considered 
method allows to solve the scattering problem in the quan- 
tum graph, to find the energy spectrum and to construct 
the wave functions. The classical analog of the presented 
vertex amplitudes method is the Kirchhoff's nodal poten- 
tials method in the theory of electrical circuits. 

Let us consider the one-dimensional scattering prob- 
lem for a quantum graph with the potentials that are 
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Fig. 1. The one-dimensional scattering on the potential well 
(barrier) as a problem on the graph. Wavy line represents the 
potential well (barrier), and ^2 are the values of the elec- 
tron wave function at the graph's vertices, 71, 72 and 73 are 
the graph's edges. The coefficients (0^,6;), i — 1,2,3 are the 
coefficients of the wave function: 4>i(x) 



placed on the graph's edges. The quantum graph is a met- 
ric graph with the self-adjoint Hamilton operator defined 
on it [35] ■ 

We assume that the quasi-one-dimensional dynamics 
takes place. The quasi-one-dimensional dynamics may be 
realized in the lowest sub-band of a very narrow quantum 
wire. The wire width in the real experiments can not be in- 
finitely narrow. But, for the quantum wire under low tem- 
perature the electron dynamics is quasi-one-dimensional 
because the higher transverse levels can not be excited. 

The considered quantum graph consists of a compact 
part connected with the reservoirs of the charge carries 
by the semi-infinite edges. We denote these edges as the 
in, out-edges. In this edges the asymptotic conditions for 
the electron wave functions with respect to the compact 
part of the graph are realized. 

In the general case, topology of the graph may be fairly 
sophisticated. In this connection the graph considered here 
is the two-dimensional array with the real potentials with 
finite support placed on the edges (see fig. [l])- 
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Fig. 2. General form of a quantum graph considered in the 
paper. Coefficients (flm, (oout, bout) correspond to the 

in,out-wa,ve functions (see equations |T] [21). Wavy lines rep- 
resent the potentials on the edges. Dashed lines show that the 
node with only two intersected edges equivalent to the one edge 
without nodes. 



are: 



The wave functions of the electron in the in, out-edges 
■>pin = ain exp(ifcx) -I- bi„ exp(-ifcx), (1) 
■>Pout = a-out exp(ifcx) -t- bout exp(-ifca;), (2) 
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where k is the electron wave-number. As it follows from ([TJ 
and ([2]) the elastic scattering takes place. This assumption 
is not necessary and does not influence on the generality of 
the proposed method, but facilitates further mathematics. 
The physical meaning of the coefficients a™, hm, Uout, bout 
depends on the considered problem (see tab. [J). 



Table 1. Coefficients of the wave functions for difi^erent "scat- 
tering problems" ("sc.pr."). 



"Left sc. pr." 


First Jost's solution 


ai„{k) = 1 


ain{k) = Mil 


bin{k) = r(k) 


bi„{k) = M21 


aout(k) = t{k) 


aout{k) = 1, 


bout{k) = 


bout{k) = 


Second Jost's solution 


"Right sc. pr." 


ain{k) — 


Clm(fc) = 


6in(fc) = 1 


h„{k) = t{k) 


aout{k) = -Mi2 


aoutik) = r{k) 


bout{k) = Mil 


bout{k) = 1 



In table [T]M represents the transfer matrix |551I57] and 
r{k), t{k) are the reflection and transmission amplitudes 
respectively [38ll39] . The explicit form of M is discussed 
below. 

In each edges 7„ of length l^^ the proper coordinates 
are used ^-y„ £ [0, In the in, out-edges coordinates 
are defined in a different way: G (— oo,0), ^o«t G 
(0,cx)). These coordinates are just a natural parameter 
in differential geometry |40) . The electron wave function 

in the graph is represented by a set of components 
'4'in = 1)2,..., M, where M is the number of edges 

in graph. Each element ) of the set is governed by 

the Schrodinger equation: 



^7,> V'l^ (^7,. ) = ^7,> (^7^ ) ' 
I d 



d 



2 (i^-y,^ \''"'7„ c'?7„ 
i?(i/^J ^{V'^^ : V7„ eCo°°[7„]}, 



(3) 



1. the standard continuity conditions at vertices u„: 

V'7l(?^J ••• =^7iv(Ci.J (4) 

where is a number of edges intersected at vertex Vn , 
^i,^ is the value of local coordinate in vertex 
and !f„ = if(C^,„); 

2. the restriction on the sum of derivatives of the wave 
functions at vertex u„: 



N 



1 dip. 



= 0, 



(5) 



where 



-1, for the edges terminating at vertex u„ 
I, for the edges emmanating from vertex Vn 



and the index j numerates the edges which are incom- 
ing at vertex or outgoing from it. 

Let us note that the boundary conditions ([5]) are similar 
to the well known Kirchhoff rules |2B,43 . 

In what follows, we work with the self-adjoint operator 
Hr only and the tilde is omitted for simplicity. 

Now, we consider the edge 7„ with an potential that is 
localized somewhere on this edge. The wave function t/;^^ 
can be determined by the following pairs of coefficients: 
(a^„ , ) and (c^„ , d^^ ) . They have the following mean- 
ing: the wave function in every edge outside the potential 
can be taken in the form aexp(ifca;)-|-6exp(— ifcx), thereby, 
the pair (a^,^ , ) are the coefficients of ij)^^ before the po- 
tential and the pair (c J are the coefficients of t/j^^ 
beyond the potential. 

It should be noted that the wave functions in the edge 
is not necessary to be a linear combination of the plane 
waves exp(±ifcx). If there any potential on the edge that 
nowhere turns into zero then the wave function is a linear 
combination of two appropriate solutions of Schrodinger 
equation. Such potentials, [e.g. V{x) = — Vbsech^7a;) were 
used for description of the aromatic molecules in the frame 
of the "graph model" (when the aromatic molecule can be 
treated as a graph) [il] . 



where rn^^ is an effective mass of the electron in the edge 
In, is the real and locally measurable potential placed 
on the edge 7„. The operators H~^^ are the symmetric op- 
erators with the defect indexes (2, 2) It is known [U 
[42] . that there exists a self-adjoint extension Hg of oper- 
ator 

M 
i—n 

where M is a number of edges in the graph Q. The self- 
adjoint extension Hg can be constructed by imposing ad- 
ditional boundary conditions on the wave functions ip^^ [U 
namely: 



3 Main formalism 

Let us consider the left scattering problem. Thus, accord- 
ing to table [U the m, out-wave functions are: 

4>in — exp{ikx) + r{k) exp{—ikx), 

^out = exp(ifca;). ' 

An application of the conditions (0 leads to the 
system of linear algebraic equations for the vertex ampli- 
tudes * = (*„ = *i, ^-2, . . . , "^out = *n)^. This system 
contains the coefficients of functions V'7„ • Hence to solve 
the obtained system for 'S' it is needed to express a^^ , 6^,^ , 
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Cy^ and in terms of (^I^i, ^'2, • ■ • , ^'n)- To do this we 
use the following obvious relations: 



(7) 

The pairs of coefficients ia-y„,bj^) and (c-y^,d-y^) are re- 
lated by the transfer-matrix M [551157] : 



The elements of the M-matrix are 



(8) 



(9) 



^7n (''7n ) is the transmission (reflection) amplitude for the 
potential placed on the edge 7„. Using relations (O, ([5]) 
and condition (O one get: 



r^i'-^") = f 1 , where 



C -C* \ -C* 1 / 

^ ^ / (10) 

C = exp(-iA:C.„+jM(]"^ - exp(ifce„„+ JM^^-^ 
From ©, (0) it follows 

r(2,7„) = [M(^„)]-ir(i'^"). (11) 

Finally, ©, ^ and ([10]), (HI]) lead to the system of linear 
equation for unknown vector yp. Finding the ^ solves the 
transport problem for graph because the wave functions 
^in, V'out satisfy the following conditions at the in, out- 
vertices: 

l + r(fc) = *„„ t{k)^^out- (12) 

The wave function on the edge can be constructed us- 
ing ([7]), (ITUl) . If there is no potential on the edge then 
r(i'Ti) and r(2.72) are identical and: 



_ -r(i.7,.) 



or 



2 sinfcL 



[exp(-ifc^J^'„ - ^'„+i] , 



2sinfcZ. 



[*„+!- exp(ifc/^J*„] 



7n 



The last gives 

V'7 



sin fc/. 



[^'n+i sinfc^^^ — sin(fc^^^ — kl^^)] . 



771 



(13) 

The obtained form of the wave function ([T3]) still correct 
for a general graph without potential on the edges and is 
used in a number of works [^[^[15] . 



The wave function (|13p creates the current j^^ in the 
edge 7„: 



= fc|*n||*n+i| sin(arg*„ - arg*„+i) cscfc^„, 

1 (14) 

CSC /c/o 



sin kl 



7n 



As it can be seen from ([T^ the current on edge is de- 
termined by the relation between arguments of the wave 
function at the appropriate vertices. This situation is 
similar to the Josephson current |46j . and, more general, 
to the interference of the currents in the graph. It follows 
from ((Ti]) that current turns into zero when one of the ab- 
solute values is zero or if arg — arg ^,1+1 — t: n, n £ Z . 

Below we demonstrate the application of the proposed 
method on some well known quantum graphs. 



4 Scattering by the quantum ring 



In this section we consider the well known example of 
quantum ring [5^IM1H7] (see eq. This model plays 

the same role in quantum transport as the harmonic oscil- 
lator in quantum mechanics. There are plenty of articles 
examined the transport properties of the quantum ring. 
The quantum ring is widely used (or proposed to be used) 
in different fields, namely spintronics [I51HU1I5D] . near-field 
optics |5T] , optics of nanostructures [55] , Aharonov-Bohm 
interferometry [531154] etc. 




Fig. 3. The scheme of scattering by the quantum ring. The 
lengths of the upper and lower arms are different, namely 
and Z-y2 respectively. The form of the arm can significantly dif- 
fer from the semicircle. The main demand to the curve repre- 
sented the arm is the smoothness. This condition provides the 
existence of the tangent to this curve. 
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An application of the condition (O at each vertex yields: 
at vertex vi: 

2 



n=l 



at vertex V2'. 

2 



(15) 



(fc) + E{[rn''^exp(ifc^„; 



n=l 



-r^''^^^exp(-ifc^j' 



^-1 



~r^2''''^exp(-ifc/^J ^2} =0 



Applying the continuity conditions (U) at vertices ^E*! and 
^'2 one obtains: 



t(fc) = *2, l + r(fc) = *i 



(16) 



Equations (fT5|) and (fTC)) represent the system of linear 
equations for the unknown vector of vertex amplitudes 

* = (*l,*2)^ 



= F, r = I + 



(17) 



where 



2 



22 ' 



(r0i2 = Eri2'^' 

n=l 
2 

(rt)2i - E ^ii"'"' exp(ifc^J + r^'i'^") exp(-iA:Z^J, 

■11=1 

2 

(rt)22 = E ^12''"' exp(ifc^J + r^'s'^") exp(-ifcZ^J, 



1 
1 



The solution of (H?]) is: 

* = (i + rt)-^F, 



*i A ^ / r22 
*2 / detr I -r2i 



(18) 



Finally, according to (jl6p . one obtains for transmission 
and reflection amplitudes: 



detr' 



detr 



(19) 



Note, that the vector 'S' can never be zero. The physical 
meaning of this fact is that simultaneous decoupling of the 
ring from the both external leads is not possible. This fact 
follows from the current conservation law. 

Equation (|17l) is the main equation of the vertex am- 
plitudes method. Once we construct the r-matrix the fur- 
ther analysis of the problem can be completely conducted 
by means of its elements. The transmission and reflection 
amplitudes, wave function of electron, probability current 
density, Green's function (not presented in this paper) are 
all expressed through the r-matrix elements. 

The mathematical base of the method can be formu- 
lated as follows: the vertex amplitudes method reduces the 
Schrodinger equation for the quantum graph to the system 
of linear algebraic equation (see eq. (IT7|)) for the unknown 
vector of vertex amplitudes * = {^'i, . . . , 



4.1 Symmetric quantum ring 

Let us begin with the most simple case corresponding to 
the symmetric ring (ring with equal arm lengths) with 
absence of any potentials on the edges. It suggests that 
M is identity matrix [55]: 

Mil = M22 = 1, M12 = M21 = 0, C = exp i-ikl) . (20) 

Since M is identity matrix then r = r^^''*'"). More- 
over, when the edges of the ring are identical the matrices 
r(i^Ti) and r(i^^2) are identical too. 

Thus, the elements of r-matrix for the ring with two 
equal arms of length I take the following form: 



11 



1 



71) p(l,7l) 
^ 21 



ri2 = 2 ri 



.(1,71) 



-r 



(1.71) 
22 



r2i 



■ 22 — 



t[\''''> exp(ifc/) - r^\'^^) exp(-ifcO 
r^2''^^^ exp(ifcZ) - r^2'^^^ exp(-ifcO 



(21) 



- 1. 



Recalling the definition of r^^'^") matrix ([TU)) gives the 
explicit form of the differences in (j2ip : 

"^11 "^21 ^ cot klry_^ J 



^(1,71) 



.(l:7i) 



— —icsc kl 



7n' 



r^'J^''*^^-' exp(ifcZ.y^) — r2i'''^'' exp(— ifc/^^) = icscfcZ^^, 



^(1,71) 



exp(ifc^^) — r22''^'' 6xp(— ifc/^^) = — icotfcZ-, 



(22) 

Therefore, the r-matrix for quantum ring with equal arms 
of length I has the form: 



l + 2icotfc/ — 2icscfc/ 
2icscfcZ — 1 — 2icotfc/ 



(23) 
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Using (|T^ and ([25]) we obtain the following useful for anal- 
ysis expressions for the transmission and reflection ampli- 
tudes: 

m = T-TnT^^7—7Tn^ (24) 



r{k) = - 



5 sin(fcZ) -I- 4i cos(fc/) ' 
3 



5 + 4i cot{kl) ' 

As it can be seen from (0^ . t{k) never turns into zero 
It is easy to see that values fc„ = ttu/I, n = 1,2, 
correspond to the zeroes of reflection amplitude. 



(25) 



there are two ways to investigate the transport through 
the ring. In short, we could fix the arms lengths vary- 
ing only the wave number and vice versa. In both cases 
condition (j29p gives the narrow (even of zero width) res- 
onances of t{k) while condition ([50)1 corresponds to the 
more smooth resonances with finite width. The first reso- 
nance will be referred to below as the first type resonance 
{FTR), the second one - as the second type resonance 
(STR). 

Let us first consider the energy dependence of the trans- 
mission amplitude (see fig. |4]). Transmission coefficient 



4.2 Asymmetric quantum ring 

The situation dramatically changes if one takes the ring 
with different arm lengths (asymmetric ring). In this case 
the F-matrix takes form: 



Fii = 1 -I- i cot fcLj -I- i cot kl. 



'72 ; 



-icscfcL 



icscfcL 



T21 — i CSC kl-y-^ -\- i CSC kl^ 



F22 = — 1 — icot kl 



71 



'■72 ) 



i cot kl 



(26) 



72- 



As before, the transmission and reflection amplitudes are 
defined by (|T^ and it is easy to show that: 



t{k) — i(sinfcZ,yj + sin klj^) x [| sin fc/,.yj sin fc^^-l- 
-|- iam{kl^-^ + kl^^) — cos kl^-^ cos kl-^^ +1] ^7 
r(fc) = I [sin fcZ-yj sin -I- i sin(fc/,^j + kl^^)] x 
X [| sin fc/^j sin kl^^ + i sin(fc/,^^ + kl^^ ) — 



(27) 



(28) 



' cos kl-y-^^ COS /t/'-y2 



II-} 



- 1. 



In the quantum ring transport experiments the ques- 
tion of resonance transmission has a fundamental aspect. 
By resonance transmission is meant not only the loss- 
less passing of electron (transmission coefficient is equal 
to one) but also the total absence of transmission when 
transmission coefficient is equal to zero. The perfect trans- 
mission (|i(fc)| = 1) is observed for electron energies equal 
to the eigenvalues of the compact graph, in our case it is 
the closed ring with different arm lengths. 

From ([27| it follows that t{k) becomes zero if 



kl 



72 



kl-y-^ = 27rri, n = 1, 2 



fc^ - kl^^ = (2n + l)7r, n = 0, ±1, ±2 . 



(29) 
(30) 

These conditions are the same as obtained by Xia [51] . 

The problem of the quantum ring under consideration 
contains the two type of parameters: energetic parameters 
represented by the wave number k and geometric param- 
eters represented by the arm lengths Z,^^, 1^^. Therefore, 




10 12 14 16 18 20 22 



Fig. 4. Different type of the transmission resonances. The 
transport through the ring can be significantly different for 
commensurate (dashed line, l^^ — Inm, l^^ = 2nm) and 
incommensurate (solid line, l^^ = Inm, l^^ = 2.1 nm) arm 
lengths. The diamonds correspond to the condition (|29|l . while 
the black points to condition (|30|l . 



demonstrates characteristic peculiarities in the case of in- 
commensurate lengths for the values of k satisfying the 
condition (|29p. Namely, the very narrow dips. Such be- 
havior can be explained by analysis of the series expansion 
of the transmission amplitude t{k). For instance, near the 
67r 

point k ~ — — , marked in fig.|l]by the diamond, it has the 
following form: 



where 



f3Ak 

Ak + n' 



Ak 



^71 ^~ ^72 



fc-6.08, (31) 



P = 2i- 



-1 + 2cos 27r 



^72 ^71 



lj2 ~l" l-fi 



sin Gtt 



lj2 + ^71 



and 



-1 + cos ( 6w 



n 



I. 



^72 ~i~ ^71 



4(^2 + ^71 ) - i(^2 - l-yi ) sin ( ^^-f^ ^"'^ 



^72 ^" ^71 , 
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The parameter in the denominator of this expan- 
sion determine the width of the resonance. One can eas- 
ily obtain that for commensurate lengths (Z^^ — 1 unij 

= 2 nm) this term is equal to zero and amplitude t{k) 
is equal to one. In opposite case of the incommensurate 
lengths (namely,Z^2 — 1 nm, Z-y^ — 2.1 nm) the parame- 
ter f2 is nonzero and in fig. |4] we see the narrow reso- 
nance. Figure|3]illustrates that the resonances width grows 
with increasing of k. In fact, the resonance width changes 
periodically upon varying the k (see fig. [5]). The reduc- 
tion of the transmission coefficient, similar to presented 
in fig. [SI was revealed by Wu and Mahler [55] in studies 
of Aharonov-Bohm effect in metals and semiconductors. 
The STR position are marked by the black points in fig. H] 



T 




Fig. 5. The dependence of the FTR width on the electron 
wave number k, l^^ = Inm, Z^j = 1.1 nm. The black squares 
show the position of the FTR, black points correspond to the 
STR. The periodical variation of the parameter \Q\ in (|31|) is 
illustrated in the inset. 



and One can see that they are significantly smoother 
compared with the FTR. There width is always finite for 
a finite values of k. 



Now, let us consider the "geometric" resonances con- 
ditioned by the dependence of the transmission amplitude 
on the ring lengths. The wave number of electron is fixed 
in this case. As it was noted above there are two different 
type of resonances (see fig. The forms of resonances 
resemble the ones considered above (see fig. S]). Never- 
theless, there is important distinction between these two 
cases. Expanding the numerator and denominator of t{k) 
defined in (P7)) in powers of kl^^ up to the first-order terms 
near two points klj.^ = 27r — kl^-^ and klj.-^ = tt + kl-^-^ we 



T 




Fig. 6. The transmission coefficient versus second arm length. 
The general view of the transmission coefficient resembles the 
one presented in Fig. U One can see the characteristic form of 
the FTR. Their positions are marked by the squares. Depend- 
ing on the upper arm length the width of these resonances 
is changing, l^-^ = 0.22 nm (solid line), Z^^ = 0.87rnm (dashed 
line). Dashed-dotted line (Z-y^ = 7r/2nm) presents the situation 
when the lengths of the arms are equal. When this condition 
takes place the transmission coefficient turns into zero (black 
diamond in the a;-axis). The electron wave number is fixed to 
be = 1 nm~^. 



obtain: 

ri\ i Ay COS. X 

Kk) ~ A (1 ■ ^ — 7-\ 1 ■ 2 ' Ay = y + x- 27r; 

Ay (i sm2a; + ij — 5 sm x 

t{k)^ 

lAy cos X 

Ay (I sin 2a; — icos2x) -I- (| cos^ ^ ~ \ cos 2a; — isin2a;) ' 
Ay ~ y — X + TT, 

(32) 

where x = kl^-^ , y = klj^ ■ The main difference of two 
expansions (j32p consists in behavior at the vicinity of the 
points X — TTTi, n — 1,2,.... It can be obtained directly 
from ((32)) that the limits Ay and x — ^ vrn do not 
commutate in the first case: 

lim lim {t{k),r{k)} = {0,1}, 

(33) 

lim lim {t(fc),r(fc)} = {(-1)",0}, 

Ay—^O a;— >7rn 

and commutate in the second. This non-commutativity is 
even more sharply manifested for coefficients of the wave 
function (ai, 61), (02, 62) 

lim lim ai = const, 

X— ^7rn Ay—^O 

(34) 

lim lim ai — ioo. 

Ay^Q x^nn 

Bulgakov et. al. [32] associate such a behavior of the trans- 
mission amplitude and the wave function in the ring with 



8 



A.F. Klinskikh et at: On the electron scattering in one-dimensional complexes: the vertex amplitudes method 



the existence of the bounded states in continuum. In de- 
tails it is discussed below. 

The noncommutative properties of the wave function 
coefficients and transmission amplitude can be understood 
using the F- matrix (l26l) . On the one hand, if kl^^ — ^^rn — 
kl-y^ then Ft = 0, on the other hand kl-y-^ — nn are the 
points of singularity of all elements of Tt . This fact shed 
light to non-commutativity of the two aforementioned lim- 
its dMD. 



If one takes A 



A, 



While kL 



27rn — fcL, make the Tt be zero ma- 



trix, the values kl^^ = nn -j- kl^^ keep diagonal elements 
nonzero. The diagonal elements in this case, ±2icotfcZ-yj, 
determine the phase of the reflection amplitude, namely 
(p = arg[r(fc) -I- 1]: 



tan(p = — 2cotfcZ 



71 • 



(35) 



It can be useful to note that at the same time ip is the 
phase of the wave function at vertex vl according to dH 



4.3 Asymmetric Aharonov-Bohm ring 

The problem becomes more complicated if at the center of 
the ring there is an infinitely thin solenoid carrying finite 
magnetic flux. At the same time it is a standard situation 
in the quantum transport [5^ . 

Let us consider the asymmetric quantum ring with dif- 
ferent arms of length and l^^ penetrated by the mag- 
netic flux From the mathematical point of view it is 
picewise-smooth contour. Thus, it can be parameterized 
by the two natural parameters, namely, and ^^j. We 
choose the "solenoid" in such a way that the vector po- 
tential A has only tangent to the contour component (see 

fig. [ZD. 




Fig. 7. The quantum loop with arms of arbitrary form cou- 
pled to the two leads. The loop is penetrated by the Aharonov- 
Bohm flux. One can see that directions of tangential compo- 
nent of the vector potential A and vector dl^.^ are opposite on 
the upper and lower arms. 



The Stokes' theorem gives 

<j> Ad\ = ^. 

L 

According to our choice of the vector potential 
j Ad\ = Ai_^J^,+ Ai^J^^. 

L 



A^ 



- A, then 



Thus, the magnetic flux modifies the phase of the wave 
function in the ring while the in, out-wave functions are 
still the same: 



ipin = exp(ifcx) + r{k) exp(— ifcx), 
V'otit = t{k) exp{ikx), 
V'71 = a^i exp(ifc~a;) + 6^^ exp(— ifc+x), 
Tp^2 — 0-^2 exp(ifc+a;) + exp(— ifc^a;). 



(36) 



where k~ = k - a, k+ ^ k + a, a ^ <f>/((f)oi), $ is the 
magnetic flux through the loop section area, $0 2TThc/e 
and L = l.y-^ -)- . 

Applying the developed in previous sections method, 
one obtains the following expressions for the elements of 
matrices T'^^'^^\ r'^in2) 



r(i.7i) 
^ 11 



T.(l:7l) 
^ 12 



T.(l:7l) 
^ 21 



.(1,71) 



T.(l,72) 

^ 11 



-.(1.72) 



7.(1.72) 

^ 21 



.(1,72) 





exp(— i/c 


''71 ) 


expi 


— ifc+^J — 


exp(ifc~^j) 




-I 




expi 


—ik+lj^) — 


exp(ifc^^j) ' 




— exp(ifc 


Z-yj ) 


expi 


—ik+l-y-^) — 


exp(ifc^^j) ' 




1 




expi 


—ik+lj^) — 


exp(ifc~^j) ' 




exp(— ifc 


'72 ) 


expi 


— ifc 1^2) — 


exp(ifc+^2) ' 




-I 




expi 


—ifc 1^2) — 


exp(ifc+^2) ' 




— exp(ifc 


'72/ 


expi 


—ifc ^72) ^ 


cxp(ifc+^2) ' 




1 




expi 


—ifc ^72) ^ 


exp(ifc+^2) ' 



(37) 



and r 

Fii = 1 + i cot fcZ-yj 4- i cot kl^^ , 

= — iexp(ia^j) cscfc^j — iexp(— iaZ^j) esc fc/721 
T21 — i exp(— ia^j ) esc fc^^ -I- i exp(ia/72 ) <2sc kl-y^ , 
T22 = — 1 — i cot kly-^ — i cot kl-^ 



'72- 



(38) 

Note that the same matrix was used by Texier and 
Biittiker in studies of the quantum ring with the Aharonov- 
Bohm flux BHl. 
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From (1551) we obtain the transmission amplitude 
for the ring with different arms , l-y^ penetrated by the 
magnetic flux (P: 



t{k) = 2i 



. exp(iQ!/^j ) sin kl^ 



sin fc/^j sin kl.y^ det F 



exp(-ia^)sinfc/^ ^^^^ 



Comparing this with (|27p. we understand that the pres- 
ence of the new parameter a modifies the conditions of the 
transmission amplitude zeros (1291) . One can see from (j38l) 
that a keeps the same points of F-matrix singularity but 
changes condition Straightforward calculation shows 
27r then 



0.8 
0.6 
0.4 
0.2 




10 12 



-k,mn 



that if kl'y2 ~1~ kl~^-^ 



121 = CSC fc/^j exp {~ial 



'72 / 



1 — exp i 



(40) 



The condition = gives $/$o = 27rn or aL = 27rn, 
where n = 0, ±1, ±2, . . .. Thus, to observe the FTR the 
two simultaneous conditions is needed to take place: 

1. kl-y^ + kl-y^ = 27rn, n = 1, 2, . . . ; 

2. aL = 27rm, m = 0, ±1, ±2, . . . . 

It can be directly verified that in the case of symmetric 
ring these conditions and ones obtained by Bulgakov et. al. [32 
appears to be the same. 

The noncommutative properties of the transmission 
amplitude emerge upon calculation the limits Ak — )■ 
and Aa — >■ 0. For the ring with equal arms it was done by 
Bulgakov et. al. [32] . Figure [5] illustrates this noncommu- 
tative behavior of transmission amplitude. 

One feature needs to be recognized. When one takes 
aL = 27rm, m = 0, ±1, ±2, . . . the value of the transmis- 
sion amplitude differs from ones with aL = 0, but trans- 
mission coefficient is the same. From the other hand, the 
values of the parameter a 7^ 27rrn can lead to the change of 
the resonance type or appearance of new resonances (see 
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Fig. 8. The nonuniform nature of the transmission ampli- 
tude in the vicinity of the FTR. The arms are of the same 
length, Z^-j = 1-^2 = 1/2 nm. The conditions of the FTR 
are formulated above and coincide with the ones obtained 
by Bulgakov et.al. [32]. It is seen that the limits Ak — >■ 
and Aa — >■ do not commutate: (a) lim lim = 0, (b) 

lim lim = 1. Here Ak — k~ ko, Aa — a — ao and fco, oo 

Aa—i-2'K Ak—^2TT 

are the points of FTR. For example, these can be fco = 27r, 
ao = 27r. 



5 Conclusion 



4.4 Scattering by the cascade of quantum rings 

Let us consider the quantum transport in the cascade of 
quantum rings. It is well known that convenient method to 
work with one-dimensional quantum cascades is the trans- 
fer matrix method |36ll37j . Transfer matrix of the linear 
chain of N quantum rings and transmission amplitude of 
this cascade are calculated as: 

M*°* =Mi...Mn, 

(41) 

t{k) ^[M\fr\ 

where Mi, . . . , Mn are the transfer matrices of interme- 
diate rings, M*°* is the transfer matrix of the chain and 
t{k) is the transmission amplitude of the chain. 

Using expressions for transmission and reflection am- 
plitudes dMl), (EH), (EZl), (EHl), (133) one can construct the 
transfer matrix of appropriate quantum ring (symmetric 
or asymmetric) according to the definition 

In fig. [To] the transmission coefficient for the different 
chains of quantum rings as a function of the electron wave 
number is presented. 



In this paper, we have demonstrated the successive analyt- 
ical approach to the solution of the transport problem for 
a quantum graph. Within the framework of the developed 
approach the solution of the transport problem is equiva- 
lent to the solution of a linear system of equations for the 
vertex amplitudes All major properties, such as trans- 
mission and reflection amplitudes, wave function on the 
graph, probability current, then are expressed in terms of 
one F-matrix that determines the transport through the 
graph. 

The simplicity and convenience of the proposed method 
have been illustrated on the example of quantum ring on 
the different assumptions. The compact part of the graph 
under consideration has the two-terminal structure. How- 
ever, the possibility of the expansion of the vertex ampli- 
tudes method on the many-terminal systems is obvious. 
This will resulted only in the increasing of the size of the 
F-matrix . 

The carried out analysis revealed the two kind of res- 
onances upon transport through the ring. The most in- 
teresting is the first one. The nonuniform behavior of the 
transmission amplitude (see l33p . features of the resonance 
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k,nm 



Fig. 9. The illustration of the induced by the Aharonov-Bohm 
flux change of the resonance type and appearance of new res- 
onances. The black squares correspond to the condition (|29p . 
the black points correspond to the condition (|3Up . Dashed line 
presents the transmission coefficient for l^^ = 1 nm, l^^ = 2 nm, 
a = 0, solid line - l-y^ — Inm, l^^ ~ 2nm, a = 3. Thereby, the 
Aharonov-Bohm conditions induce the FTR even for commen- 
surate arms lengths. 



width, the specific structure of the resonance in the com- 
plex energy plane (not presented here) make the FTR very 
similar to the BIC considered by Pursey and Webber [S^ • 
In this paper wc intentionally avoided the use of term BIC 
for the quantum ring. From our point of view, the direct 
analogy with the results of von Neuman and Wigner [55] , 
Pursey and Webber |57] is not correct in the case of the 
quantum ring without any potentials. One of the central 
idea of von Neuman and Wigner was in special choice of 
the potential the electron is subjected to. This choice leads 
to the integrability of the continuum wave function. The 
relationship between the properties of the wave function 
and potential is crucial in the definition of the BIC. For 
the quantum ring we observe the emergence of the FTR 
(in both cases of rings with and without Aharonov-Bohm 
flux) in the absence of any potential. 

The transmission resonances for the quantum graphs 
and, in particular, for the quantum ring, were studied in 
a number of works [53l[58l[59ll60] . Upon studying the res- 
onance line-shape features it was established these reso- 
nances are the Fano resonances. The Fano resonances ap- 
pear when the coupling between the discrete spectrum and 
continuum takes place. Most often, the discrete spectrum 
originates from the transverse quantization. However, even 
in a more simple quasi-one-dimensional model the Fano 
resonances emerge [SS]- In this case, they come from the 
coupling of quasi-bound states of the compact part of the 
graph and continuum of the external leads. This in turn 
gives rise to the peculiarities of the the resonance struc- 
ture in the complex energy plane [58] . The classical analog 
of this phenomena are the Helmholtz resonances [HI] . The 
ring play the role of the resonator connected to the waveg- 
uides. The quasi-bound states can leak out from the ring, 
that is the coupling is done by the contacts. 

Let us notice that the theoretical framework develop- 
ing during research of quantum graph is interesting for 
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Fig. 10. This figure shows the transmission coefficient through 
the linear chain of the connected quantum rings. The line- 
shape can be rather different depending on the ratio of the 
arm lengths: (a) transmission coefficient for the three quan- 
tum ring with different, but commensurate arm lengths l-^-^ = 
1 nm, = 2 nm; (b) transmission coefficient for the three 
quantum ring with different and incommensurate arm lengths 
= Inm, — 2.1 nm. 



a field that are, at first glance, far from the considered 
problem. These are neurobiology and cellular biology. The 
first reason follows from the usability of the graph to es- 
tablish the relationships between the different parts of 
complex system. The modular constitution of biological 
systems becoming more and more convenient [62) . Thus, 
the action potentials characteristic of nerve and muscle 
cells have been reconstituted by transplanting ion channels 
and pumps from such cells into non-excitable cells |63) . 
Thereby, the properties of single modules and molecular 
connections between them are similar to electrical circuits 
(with the Kirchhoff rules) and quantum graphs (with the 
principles of quantum mechanics). 

The second reason concerned with the existence of bi- 
ological systems that are the natural graphs. For example, 
it is neural networks in brain 64,65,66 . The fundamen- 
tal properties of neural networks strongly depend on its 
topology and coupling conditions at the nodes. Note, that 
the same is typical for the quantum graphs. 
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A Scattering by the potential with finite 
support 

By way of illustration of the vertex amplitudes niethod, 
let us consider a well known problem of electron scatter- 
ing on the one-dimensional potential with finite support 
(see fig. [TT]) . Below, we derive some general expressions 
for transmission and reflection amplitudes using the ver- 
tex amplitudes approach. It must be pointed out that the 
obtained expressions are valid for the potentials with arbi- 
trary profiles (not only for the square barriers, wells etc.). 
Let us consider the left scattering problem in the infi- 



Yi 72 73 
mmimi* 

Fig. 11. Schematic representation of the scattering by the po- 
tential with finite support. Wavy line represents the potential. 



nite quasi-one-dimensional lead with real potential V{x) S 
. In fact, these requirements are more stringent than 
those actually needed for the validity of the results re- 
ported below. Following the graph ideology we treat this 
system as a graph consisting of three edges 71 , 72 , 73 and 
two vertices fi, W2. The local coordinates € (—00, 

S [C-ui ] ^-iid e [^ij , -fcx)) are defined in their 
respective edges. Let potential V{x) to be located on the 
edge 72. For simplicity, we assume that the electron's ef- 
fective mass m-y^ to be equal to the mass of the free elec- 
tron. According to the scattering theory ideology the wave 
functions in tlie edges 71, 73 are taken in the form (p]). 

According to theory of differential equations [57] , there 
are two linearly independent solutions u{x) and v{x) of the 
Schrodinger equation in the edge 72. Thus, the wave func- 
tion Tp-y^ can be chosen as a linear combination of functions 
u and V, namely: 



^'72(^72) = a u(^^J + 6^(^72 )• 



(42) 



Let us express the coefficients (a, b) of the ip-y^ in terms of 
iZ'i and iZ'2, where 'Pi and iZ'2 are the values of the electron 
wave function at the vertices vi and V2 respectively. To do 
this, we use the obvious relations: 



Solving this system one obtains: 

1 



(43) 




(44) 



Applying the hermiticity conditions ([5]) to the wave func- 
tion at vertices wi, V2 one obtain the two equations: 

ikexpiikCJ)[l-r{k)]-[au'{^l^^^) + bv'{C,^)\ =0, (45) 

- ikt{k) exp(ifcf(3)) + [au'i^if) + bv'i^if)] = 0, (46) 

where a prime denotes differentiation with respect to ap- 
propriate local coordinate. Using the continuity condi- 
tions (HI in vertices and (H^ we obtain the system of two 
linear equations for unknown vector of vertex amplitudes 
^ = (vti, \I'2)"'" that can be written in the matrix form: 

r* = F, (47) 

where 



rii 








-fifcA, 


ri2 


= u{df) 




-vidfwidf), 




r2i 






-u'idfHdf), 




r22 






-uidfwidf) 


-fifcA, 



^ ^ f 2ifcAexp(ifc^ii^) 
V 

The solution of (gZl) is: 

* = (ifcAi + rf)-^F, 

2ifcexp(ifc^i,i^)A / F 



detr 



22 

-r2i 



(48) 



The transmission and reflection amplitudes are easily cal- 
culated according to ^ 



(49) 



t{k)='f2CXp{-ik^if), 

r{k) = [*i -exp(iA:Ci^))] exp(ifce^J)). 
Substitution of dM]) into (|49]) gives: 

2ifc(i.2-M2)exphifc(^g)-e^)i^)] 



m = —I 



(v2 — '^k){pLi + ik) — —{vi — ik){^2 + ik) 

U2 V2 



where 



(50) 



(ui,2; 1^1,2) = (uidf.2);v{df.2)) , 



(Aii,2;i^i,2) 



? 5^1 2 
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The expression for the reflection amplitude can be ob- 
tained in analogous manner, namely: 

^j.^^ _ M2fi(/i2 - ifc)(t^i - ifc) - uiV2{v2 - ifc)(Mi ~ ifc) 
uiV2ii^2 - + i^') - "2^1 (j^i + ifc) (^2 - ifc) 

X exp ^2ifc^ii^^ . 

The transfer-matrix of potential V(x) could be written 
down according to 

It is significant to point out that the F-matrix carries 
information about discrete spectrum of potential V{x). 
Using the standard method [35] , or directly from the prop- 
erties of the transmission amplitude (SS] it follows: 



detr(fc) 



A(fc) 



= 0, 



(52) 



K„ — 1m\En\l'h , En is the energy of the n-th bound 
state of the electron in potential Vix). To explain it, let 
us note that the denominator of expression ([50]) is equal 
to detr/A(fc). Therefore, the solutions of (15^ correspond 
to the poles of the transmission amplitude ([50)) . 

Thus, the proposed technique allows to resolve the 
scattering problem in terms of two linearly independent 
solutions of the Schrodinger equation for potential Vix) 
and to calculate the binding energies of this potential. 



B Scattering by the system of parallel 
quantum wells 

The problem of electron scattering on the system of n 
connected in parallel quantum wells (see fig. I12[) can be 
easily solved using the vertex amplitudes method. To do 



following form: 

n 



p(l,7a) 

21 ' 



ri2 = Er 



(1,7=) 
12 



r2i = J2 ^11^^^ exp(ig,I^J + 4^1'^=) exp{-iqJ^J, 

s=l 

n 

r22 = J2 ^12^°^ exp(ig,^) + r('2'^=) exp(-i<z,^). 



s=l 



Here Qs — \/lt^~-Vs, Vg = 2mJJs/TT? and Us is the well's 
depth on the edge s. Let us transform to dimensionless 
variables by taking 1 nm as length's unit. As one can see, 
the above F-matrix coincides with the F-matrix of the 
quantum ring up to the replacement of q by fc. If all wells 
have a similar depth Uq and width /, then the expression 
for F-matrix can be significantly simplified. Using 
one obtain: 



/ ^ .nq 

/ 1 + 1— cot ql —i-r- CSC 



. nq 



.nq 



-1 — 1 — cot ql 
k 



ql ^ 



(53) 



where q — y/k'^ — Vq- The transmission and reflection am- 
plitudes can be found using equations (jl9p . Namely, for 
the transmission amplitude we get: 



i„(fc) 



2inkq 



(fc2 + 71^ q"^ ) sin ql + 2inkq cos ql 



(54) 



Y„ 




Fig. 12. The graph which corresponds to the scattering prob- 
lem on the system of n connected in parallel quantum wells. 
Wavy line represents the potential (in the considered case the 
potential is the rectangular quantum well). 



As it should be, at n = 1 this expression turns into the 
standard expression for the scattering amplitude for the 
one-dimensional single quantum well. Using (1541) we can 
analyze the dependence of the energy spectrum and trans- 
mission coefficient of such a system on the number of con- 
nected quantum wells. 

At first, let us consider the transmission probability. 
From fig. [13] it follows, for the connected in parallel quan- 
tum wells the characteristic oscillations of the transmis- 
sion coefficient take place. These oscillations are resulted 
from quantum interference effects in analogous to the above 
considered quantum rings. It is also seen that in the region 
of energy far enough from the well's bottom, k ^ -^Z— Vb, 
the period of oscillations does not depend on the number 
of wells. Indeed, the elementary calculations show that 
in this case the transmission coefficient becomes unity at 
fc — TTs/l, s = 1,2,... and minimal at fc = (2s + l)7r//, 
s = 0,l,2,.... The minimum value Tmin is approximately 
determined by the following expression: 



this one need to construct the F-matrix of the system. 
In the considered case the elements of F-matrix have the 
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Fig. 13. Transmission coefficient vs. wave number of electron 
k. Solid line corresponds to the single well, dashed line cor- 
responds to the two wells, dashed-dotted line corresponds to 
the twenty wells. All quantum wells have the same depth and 
width, namely Uo ~ — 0.5eV and / = Inm. The single quan- 
tum well with of such depth and width possesses one bound 
state with energy equal to — 0.04eV. 
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Fig. 14. The logarithm of the transmission coefficient as the 
function of the number of identical parallel quantum wells. Ad- 
dition of the new quantum well leads to the shift of the bound 
state energy position closer to a limiting value n = \/—Vo. 
Solid line corresponds to single well, dashed line corresponds 
to three wells, dashed-dotted line corresponds to ten wells. The 
single well parameters are Uo = — 0.5eV and I = Inm. 



This expression shows that with increase of the the num- 
ber of parallel quantum wells the resonances of the trans- 
mission coefficient become "narrower" and values of trans- 
mission coefficient between resonances become closer to 
zero. 

Using (IM)) we can analyze the energy spectrum of the 
system under consideration. To do this, recall that ener- 
gies of bound states correspond to the simple poles of the 
transmission amplitude .38..39) . In order to determine the 
bound states energies construct the analytic continuation 
of the transmission amplitude t{iK). The bound states are 
placed in the half-line Imfc > 0. In fig. [14] the maxima of 
the logarithm of the transmission coefficient square corre- 
spond to the bound states of the system of n connected 
in parallel quantum wells. With increase of the number 
of the quantum number the bound state energy increases 
and tends to a limiting value k = \/—Vo, fig. [T51 

It is interesting to note that if one change the number 
of connected in parallel identical quantum wells the new 
bound states in the system will not appear. 
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Fig. 15. The dependence of the bound state energy in a system 
of parallel coupled quantum wells on the number of wells n. 
The dashed horizontal line shows the limiting value k — y/—Vo- 
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